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Gravitational production of KK states 
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Gravitational particle production in the context of braneworlds is considered from a phenomeno- 
logical point of view. The production of KK modes for bulk fields is discussed and their abundances 
computed. The results have been applied to some specific fields such as gravitinos and axions, and 
their cosmologies have been outlined, exemplifying the constraints on the properties of the extra 
dimensional model in each case. 



I. INTRODUCTION 

Gravitational particle production in an ordinarily ex- 
panding Friedman-Robertson- Walker (FRW) Universe 
dominated by dust-like or radiative matter is known 
to be a very poorly efficient particle creation mecha- 
nism Indeed, if one leaves the conformal anomaly 
exception aside, despite being incessantly operative dur- 
ing the whole evolution of the Universe, the abundances 
of non-conformally coupled particles is totally negligible. 
This is due entirely to the weakness of the coupling of 
particles to gravity, combined to the steadiness of the 
evolution of the Cosmos. In order for this particular 
mechanism to be efficient one needs to look for abrupt 
transitions in the dynamics of the FRW scale factor, such 
as that believed to have taken place during inflation, and, 
more importantly, when the de Sitter phase drew to an 
end, repopulating the Universe. 

However, as it has been recently pointed out in ref. 0, 
if the four dimensional FRW Universe were to be embed- 
ded into a higher dimensional spacetime, such as in 
then even dust-like or radiation dominated FRW Uni- 
verses would be able to inject gravitationally produced 
particles into the plasma in a sizeable way. This is easily 
understood as a consequence of the possibility that the 
actual scale of - 5D - gravitational interaction be tuned 
to much lower values, thereby drastically enhancing the 
coupling with matter. Hence, for the same reason one 
expects to see signatures of (collections of) Kaluza-Klein 
(KK) gravitons in TeV colliders m the most opti- 
mistic scenarios, it is expected that expanding cosmolog- 
ical backgrounds should turn out to be phenomenologi- 
cally relevant as well. 

This new feature of brane cosmologies shows up when 
one derives the effective Friedman equation for a 4D 
brane immersed into a 5D spacetime. In this case it is 
known that the evolution of the Universe changes at high 
energies, allowing for a new term in the Friedman equa- 
tion for the Hubble parameter, which in general reads Q 
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brane, M4 ~ 2.4 • 10^^ TeV the 4D reduced Planck mass, 
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is the brane tension and M5 the effective gravity - reduced 
- scale of the five dimensional theory. The transition tem- 
perature Tt is the temperature at which the evolution of 
the Universe switches from braneworld regime to stan- 
dard one and, if the Universe is radiation dominated, it 
is given by 



Tt 



360 



1/4 



Ml 

Ma 



1/2 



(3) 



where counts the relativistic degrees of freedom tightly 
coupled with the thermal plasma. 

Therefore, if the temperature of the plasma is higher 
than the transition temperature, the Universe would ex- 
pand according to H = p/6M| . One can immediately see 
the effect of the 5D embedding: at high energies the ex- 
pansion is regulated by the 5D Planck mass rather than 
the low energy 4D coupling. 

At this point it is useful to review the basic results 
found in which provide the basis for the extensions 
presented hereafter. In the following formulas everything 
will be expressed in TeV units, unless explicitly stated. 

As already specified, for the time being the focus will 
be on fields non-conformally coupled to gravity, whose 
gravitational 'charge' is their mass m. Gravitational par- 
ticle production in a time- variable metric is in general ef- 
ficient as long as the Hubble parameter H is much larger 
than the particle mass, that is H ^ m. In this regime it 
can be shown that the number density of created particle 
is constant H 
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where p is the energy density of ordinary matter on the 



whereas for H <^ m particle creation is negligible and n 
decreases as l/a^, where a is the FRW cosmological scale 
factor. Of course this result is slightly model-dependent, 
i.e. it depends on the specific law linking expansion rate 
to time, but numerical results confirm that it is accurate 
within an order of magnitude for any relevant power law 
expansion of the scale factor, a{t) cx f^. Consequently, for 
the sake of simplicity, for H > m, the particle number 
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density will be given by eq. ^ , and particle production is 
taken to cease instantaneously as soon as H = m. Men- 
tion to the actual numerical result will be given through- 
out the paper. 

Within this approximation, the freczing-out of gravita- 
tional production corresponds to a temperature around 
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Were this species a stable one and very weakly interact- 
ing with ordinary matter, so that it would not thermalise 
or decay as the Universe cools down, then its relic density 
today would contribute with a fraction of 
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to the total energy density of the Universe, see 0, 
where h is the present Hubble parameter in units of 
100 km s^^ Mpc^^. Eq. ^ is derived assuming that 
at the temperature Tf there are about ~ 100 rela- 
tivistic degrees of freedom, which implies also a dilution 
factor of 0.04, and is valid only if the conditions 



Tf > Tf M, < 6 • 10^ ml/3 , 



Tf < M5 



A/5 > 0.2 TO, 



(7) 
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are fulfilled, as required by consistency. 

This is the result which was used in [2| to assess the 
importance of such mechanism in, e.g., the generation of 
sterile viable Dark Matter (DM), baryogenesis scenarios 
at low temperature (see also 0): or its potential danger- 
ousness, had the produced particles been equipped with 
a decay channel timed after the synthesis of the light 
elements (Big Bang Nucleosynthesis, or BBN), as for in- 
stance for the gravitino case. 



II. KALUZA-KLEIN TOWERS 

The situation hitherto described applies to brane- 
trapped fields, which therefore depend only on four co- 
ordinates, and have by definition only a zero 4D mode. 
However, especially in the gravitino case, it is likely that 
fields have access to the full dimensionality of the space- 
time, and that, as long as the extra dimension is compact, 
from a 4D observer they would be seen as a tower of KK 
states. Such a configuration was indeed studied in con- 
nection with the thermal production of KK gravitinos for 
a number of supersymmetric (SUSY) extra dimensional 
models ^. 

The aim of this section is to extend the analysis given 
in to the full tower of KK states, whereas the following 
ones will review its implications for the phenomenology 
of extra dimensions, and its constraining power applied 
to specific fields, one at a time placed in the higher di- 
mensional spacetime. A general, theoretically oriented 
discussion about gravitational particle production for 5D 
fields can be foimd in ref. M. 



Flat fifth dimension 

In this first scenario the extra dimension will be taken 
to be flat. In general the mass of each KK mode is 
given approximately by to„ = toq + n/R, where R^^ = 
2ttM^/M^ is the inverse of the size of one compact ex- 
tra dimension. This expression of course is not univer- 
sally true, since an effective potential for the 5D field 
may appear upon dimensional reduction, and therefore 
contribute to the effective mass of the KK mode. Nev- 
ertheless, in general it is expected that at high n this 
parametrisation provides a good approximation, as it will 
be clarified in the worked examples considered below. 

Having obtained an explicit formula for the mass, one 
can readily write down the abundance associated to each 
KK mode. In the limit for which the zero mode is light 
compared to n/R, which again is taken in light of the 
particular models to be dealt with, the fraction of energy 
density and the yield variable turn out to be 
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YJl ~ 1.3•10-^ln9/4M5'/^ 



(9) 
(10) 



Clearly, eq. ([9]) holds only for stable particles or for par- 
ticles whose lifetime exceeds the age of the Universe. 
YJ^ ~ n^/s, where s is the total entropy density of the 
Universe, is a conserved quantity for species which are 
decoupled from the thermal plasma and have not yet de- 
cayed. We notice that we are still assuming a Universe 
dominated by rclativistic matter, where the amount of 
particles produced by the time-varying metric is small 
when this generating mechanism is efficient. These two 
quantities are the most useful in obtaining constraints 
on the parameters of the theory, and will be extensively 
employed. 

It is immediate to rewrite the consistency condi- 
tions ([71 [5]) for this specific setup: they are 



Tf >Tt ^ n > 5 . 






Tf <Ah nMl < 5 - 10^°. 



(11) 
(12) 



What the first line means is that the results jU [T0| do 
not hold for the first handful of modes, as the quantities 
and have been deduced assuming H = p/6M^, 
while for small n we would find Tf < Tf. For these light 
modes, the gravitational production stops at the transi- 
tion temperature Tt or, if very light, during the standard 
expansion, thereby rendering their abundances negligi- 
ble. However, it should be stressed that is the case only 
if the zero mode mass is small compared to the mass 
gap ~ 1/R between the modes. Indeed, in that situa- 
tion a heavy zero mode would be copiously produced, 
and so would be each KK mode from the very bottom 
of the tower. Anyway, since such a configuration would 
not significantly alter the results presented below, this 
possibility is not discussed further. 

Moreover, we stress again that our discussion is reli- 
able for large n's, while the zero and the first few modes 
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would require a separated treatment, which depends on 
the exact model under investigation. For example, even 
a massless zero mode may be gravitationally produced, 
because in a more general setup the 5D metric 

ds^ = e2'^(^)a2(7?) [dri^ - dx^ - dxj - dxj] - dz^ (13) 

is not conformally flat (for more details, sec the discussion 
in sec. 8 of ref. [9|). In a radiation dominated Universe, 
that is the case of interest in the present work, there 
is an accidental recovery of conformal invariance, hence, 
no gravitational production of the zero mode via time 
variation of the metric (as long as the conformal anomaly 
is not considered, see appendix). 

In general each mode will have its own interaction 
strength(s) (if any) with Standard Model (SM) fields, 
and the history of each mode would depend on the way 
they decay. The same statement applies to the impact 
on cosmological observables such as DM energy density, 
light elements abundances, and so on. The main discrim- 
inator among these different effects is the lifetime of the 
particle, which, once the zero mode's decay properties 
have been specified, depends almost solely on the mass 
of each higher mode. It is thus convenient to split the 
KK tower into various bands to be taken into account 
separately. 

Once the interesting lower and upper limits have been 
identified, one is in a position to scrutiny the number 
of modes in that given range, and investigate their to- 
tal impact on a particular observable. Moreover, it is 
reasonable to expect that only large number of states for 
each band could provide significantly different constraints 
compared to an ordinary 4D particle with mass within 
the given band. This last consideration, combined with 
the fact that the zero mode is supposed to be light, im- 
plies that the highest KK mode N in the band will satisfy 
N ^ Nq, with Nq the lightest mode in the same band. 
Then it is straightforward to obtain the overall density 
fraction and yield variable for each band by just sum- 
ming over the modes and then discarding the (smaller) 
contribution of the lowest state, the result of which being 

nxh^ = Y.^xh^ - 7-10-9iiVi^/*Mf/', (14) 

n 

n 

We note in passing that the actual summation ap- 
proaches the numerical results and -obtained 
by integration, for a number of modes N > 20, beyond 
which the precision is better than 10%. 



Warped fifth dimension 

If the extra dimension is warped, the structure of the 
KK tower turns out to be different. In the case of two 
branes, typically one finds that the mass of the KK mode 



n is 

TO„ = mo + kxne~'^^^ , (16) 

where fc = (1 - cxp(— 27rfci?))Af|/Af| is the ^^6*5 curva- 
ture, Xn is the n-th root of the first order Bessel function 
Ji and ttR is the size of the orbifold. The mass splitting 
is Am ~ 3fc exp(— 7ri?fc), and for later convenience we de- 
fine the function F{kR) = (1 — cxp(— 27rfci?))/ exp(7rfci?). 

The Friedman equation is basically the same of the case 
with a fiat extra dimension and we can still use eq. ([T]). 
It is thus fairly straightforward to repeat the steps un- 
dertaken for the fiat case, with results 

r!j/l2 ~ 3 • 10-^1 7.13/4 J7l3/4^,^15/2 ^ (17) 

~ 2.10-^2^^/4^9/4^/5'/', (18) 

in place of eqs. ([OKTO]), and 

n^h^ ^ 7-10-92^"/4^i3/4^ff/2^ (19) 
YJl ~ 8-lO-''N^^/*F^^*M^^\ (20) 

which now replace eqs. (fT4l [T5 |) . 

It is now easy to see why the scenario is more or less 
unchanged in this picture, unless one employs very small 
values for F{kR). Indeed this function is always smaller 
than ~ 0.4, for which value the abundances and consis- 
tency constraints turn out to be very similar to those 
previously obtained for a fiat extra dimension. If one in- 
stead chooses to work with much smaller F{kR), such as 
for kR =11, then the mass splitting becomes extremely 
tiny, unless the five dimensional Planck mass is pushed 
all the way up to the 4D M4. Therefore, the mass gap 
drops by a factor of 2/F, and the abundances of pro- 
duced particles increase by a huge factor 1/F (this can be 
seen explicitly once the expressions for N -see section lllll 
is plugged in (jl9|) and ([20]) ). The compensation in A/5 
is proportional to /'^4/2i again for kR ~ 11, is 

around 700, although in highly warped models it is cus- 
tomary to safely (as far as this mechanism is concerned) 
take A/5 ~ A/4. 

We would further like to stress that all our results are 
based on the validity of eq. (jT]). Actually, eq. ^ is de- 
duced from the 5D Einstein equations in the case i) there 
is no fiux along the fifth dimension and ii) the bulk en- 
ergy density is constant, which is the case for the stan- 
dard scenarios where particles arc confined to the brane 
and the energy density of the bulk is given by the 5D cos- 
mological constant. Here these assumptions cannot be 
rigorously true, but departures are usually small: even 
if there is exchange of energy between the bulk and the 
brane, in the situation of quasi equilibrium there is no 
net fiux, while the energy density of the bulk is at least 
almost constant, because dominated by the 5D cosmo- 
logical constant A5 ~ A/|. 

Along these lines one can now transfer to the warped 
solution the results which will be presented below for the 
flat scenario. 
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III. APPLICATIONS 

As already mentioned, the constraints that can be ex- 
tracted from the cosmology of KK particles gravitation- 
ally produced depend crucially on the detailed structure 
of the interaction such species have with ordinary fields. 
It is consequently mindful to analyse few cases explicitly. 
The first possibility to be taken into account is that of 
the gravitino field, simply because, once SUSY has been 
turned on, it is expected to necessarily possess a tower of 
KK modes, being part of the same multiplet the graviton 
is, therefore having access to the same spacetime dimen- 
sionality. One other possibility is gravitationally created 
axions, and will be discussed thereafter. 

KK Gravitinos 

The gravitino is the fermionic superpartner of the 
graviton in SuperGravity (SUGRA) theories, which ac- 
quires a mass via the SUSY version of the Higgs mech- 
anism. Its properties are uniquely defined by SUGRA, 
the only exception being the value of its mass, which de- 
pends on the particular mechanism that breaks SUGRA. 
Since this field is partnered with the graviton, once one 
expands the spacetime to five dimensions, it is natural 
to expect the gravitino to live in the fifth dimension as 
well, which in turn means that from a 4D observer it will 
appear as a tower of KK excitations. 

In the case of the 4D gravitino, its lifetime is a func- 
tion of its mass alone (as long as unstable gravitinos are 
concerned), being the coupling universally determined 
by the gravitational one, which also provides departure 
from thermal equilibrium throughout the expansion of 
the Cosmos. Hence, constraints on gravitino abundances 
in the early Universe are typically given for a certain mass 
range. 

Light gravitinos, with mass TO3/2 below 100 GeV, are 
usually supposed to be the Lightest Supersymmetric Par- 
ticles (LSP's) and therefore are surviving forever due to 
7?.-parity conservation. Anyhow, even if this symmetry is 
badly broken, some of the lighter modes will still decay 
with a characteristic lifetime longer than the Universe's 
age. It is mandatory to ask then that such gravitinos do 
not overwhelm the Universe, in fact, they must not ex- 
ceed the present day's DM energy density QuMh^ r:; 0.12. 

If the gravitino mass is in the range 100 GeV - 30 TeV, 
its decay can spoil the successful predictions of BBN: 
for small masses it catalyses the photodissociation of the 
light elements, while for large masses it injects pions into 
the plasma, altering the neutron to proton ratio at the 
onset of the BBN. Here the limits on gravitino abundance 
depend on the open decay channels and their branching 
ratios, but they usually are the most stringent ones. If 
the main decay mode is hadronic (Bh ~ 1), the con- 
straints on the gravitino abundance 13/2 are roughly at 
the level 10^^^ — 10^^^, depending on the exact value of 
7713/2- On the other hand, if the main decay mode is pho- 



ton + neutralino, because the gravitino is lighter than the 
lightest coloured superparticle, is about 10^'^ and the 
constraints on ¥3/2 are weaker, of order 10~^^ — 10"^^. 
For more details, see ref. [Tsj . 

The third option, although not easy to arrange in 
model building with soft SUSY breaking, is that the grav- 
itino be much heavier than a few TeV. That implies a 
lifetime shorter than about 10~^ s and no late time en- 
tropy release. However, if 7?,-parity is conserved, then 
for each gravitino that decays there will be an LSP at 
the end of its decay chain, which is then defined to be 
absolutely stable. Its abundance will refiect the initial 
gravitino one, up to 7773/2 ^ 10^ TeV (heavier gravitinos 
decay at T > 10 GeV, which permits thermalisation of 
the generated LSP's), and must be kept under control in 
order to avoid confiict with the observed DM density. In 
this case, the bound on Y3/2 is 

y3.<6.io-(i2^). m, 

V T^LSP J 

where ttilsp is the LSP mass, and is independent of the 
exact gravitino mass. 

So far for the ordinary 4D gravitino, whose abundance 
produced by gravitational interaction was investigated 
in 0. However, the picture can be dramatically mutated 
by the come into play of the KK states. In principle 
these states will have interaction properties other than 
the zero mode ones, being inherited from those of the 5D 
SUGRA model, most relevantly from the mechanism of 
SUSY breaking. The mass matrix for the KK states is 
in general non diagonal, leading to mass shifts possibly 
different state by state. For the sake of clarity though, 
and since the results outlined below are easily adapted 
to different models, each mode will be given the zero 
mode coupling, and the masses will be simply taken to 
be rrin = 7773/2 + n/R. With the details of KK gravitinos 
interactions at hand, the bounds on their abundances are 
readily computed. 

The first band corresponds to the lightest KK graviti- 
nos, with a mass smaller than 100 GeV. Here we assume 
that the zero mode is much lighter than 100 GeV and 
thus 777„ fa n/R. Under this hypothesis, the number of 
lighter states is TV ~ 9- lO^^ /M|. The total abundance of 
light gravitinos today, and its relative constraint, would 
then be 

^3/2/72 _ 8-10^^M~^^^'^ <0.12, 
^ M5 > 3-10^ TeV . (22) 

Notice that this limit assumes stability for all the modes 
included in the summation. This may not be true if there 
are allowed transitions between different KK modes, not 
taken into account in this simplified exemplifying analy- 
sis. This is expected to be a good approximation for the 
light modes since these decays have typically narrower 
widths than the usual zeroth order gravitino [3], and 
are therefore expected to be more long lived. 
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The second interesting band encompasses KK modes 
with masses 100 GeV < m„ < 30 TeV. Here the number 
of KK excitations is TV ~ 3 • 10"^^ /M| and provides the 
foUowing hmit 

8 • lO^Hl-''^' < 10-12 , 
Ms > 1.4-10*^ TeV. (23) 

This bound is actually relatively conservative, as within 
the band we have included also the more constraining 
1 TeV gravitinos, which arc typically (depending on 
their hadronic decay branching ratio) required to have 
Y3/2 ^ 10~i^. Indeed, if the main decay channels were 
hadronic one could pull the limit on down to 10"^'^ 
which immediately pushes A/5 up to 3-10* TeV or higher. 
In any case, these limits are only illustrative, as there 
are more uncertainties affecting the precise value of ¥3/2 7 
such as the more realistic numerical results for the num- 
ber density of particles (which in the braneworld expan- 
sion regime is close to a factor of 10 more), and the pos- 
sible presence of more thermalised relativistic degrees of 
freedom at that epoch. 

Furthermore, here again there could be decays which 
violate the KK-number such as KK'' KK"" + X, 
where m < n and X is a Standard Model (SM) par- 
ticle. However their amplitudes is often negligible [l3| 
with respect to the standard decay channels (since also a 
transition KK LSP + X violates KK- number). This 
argument applies to the forthcoming paragraph as well. 

The third and last band refers to the highest KK states. 
The upper bound can be taken to be of order 10^ TeV, 
leading to TV ~ 9 • 10^''/Af|, and to a total Y3/2 which 
is about 10^1 times the one found in ((23|) . that in turn 
translates into a higher lower bound for M5, that is, 

Y,/2 ~ 1.6 • lO^Hl-''^' < 10-12 , 
=> M5 > 2 • 10^" TeV . (24) 

These limits are valid for A/5 < 3 • 10^^ TeV for in this 
case we have much more than 5 KK modes in the band. 

So far it has always been assumed that the Universe, 
in the course of its evolution, has had available high 
enough temperatures, and fast enough expansion rates, 
for efficient gravitational particle production. This point 
is directly connected to the epoch which succeeded the 
inflationary expansion regime. This exponentially fast 
growing regime, besides providing another extremely well 
functioning source for gravitational creation of fields (al- 
though immediately exponentially diluted) is believed to 
have been followed by an abrupt phase transition which 
repopulated and thermalised the empty Universe, where 
again gravitational interaction could have played a major 
role. The temperature at the onset of the equilibrium era 
is usually calculated assuming instantaneous and efficient 
conversion of the energy stored in the field which drove 
the expansion, to that of the newly born SM (and its 
extensions) particles. Although not precise, this method 



gives an idea of the scales involved, and is readily worked 
out. 

As a brief example, one can consider chaotic inflation 
scenarios on the brane. The reason for this choice, be- 
yond its simplicity and its being favoured by observa- 
tions, is that the completion of the reheating process 
happens in the brane expansion regime for most values 
of A/5. The potential for the scalar inflaton field is given 
by y = m|02/2 and, as one learns from [Tsl . [T6| . suc- 
cessful generation of density perturbations in the brane 
regime requires ~ 5 • 10"^ A/5, unlike the usual 4D 
models where the inflaton mass turns out to be close 
to the 4D gravity scale, as expected. If the scalar field 
interacts with fermions, its typical decay width can be 
parametrised by F ~ 94,^^, which gives a regime re- 
heating temperature (obtained from equating F to H) 

Tr ~ Omg'/^Ah >Tt ^ A/5 < 2 • lO^^gi/^ ^ 

where the last step is required for self consistency. The 
highest mode available is therefore ~ 5 • 10'^^ g^/M^, 
obtained by requiring that Hn ~ mAr, and the overclo- 
sure constraints reads 

Y3/2 9.lo"4^/x-'<lo-'^ 

A/5 > 10i2gi3/8 TeV . (25) 

Notice that this limit makes sense only if the modes pro- 
duced do not weigh more than 1713/2 — 10^ TeV, as is 
always the case for < 0.1, for the same reasons ex- 
plained above (above g^ ~ 0.1 eq. ((24)1 applies). Lastly, 
the constraint (j23p coming from BBN is modified into 

M5 > 10"gi^/^ TeV , (26) 

using the hadronic bound, valid below ~ 2 • 10-"^, and 
the lowest band limit ([22ll becomes, up to ~ 7 • 10—*, 

A/5 > 6 • lOi^gi^/" TeV . (27) 

It is then possible to combine all the constraints in one 
single plot with the 5D Planck mass against the coupling 
constant g^. By looking at the figure [T] one sees how 
the lower bounds on A/5 are significantly weakened once 
a relatively low reheating temperature is taken into ac- 
count, still leaving interesting constraints on the param- 
eter(s) of the extra dimensions. When the inflaton cou- 
pling constant grows, the limits (HH), and (|24|) are 
approached as expected. Other inflationary models are 
briefly mentioned in the appendix, where it is explained 
why chaotic inflation is probably the only relevant one to 
study in this context. 

There is one more subtlety to be discussed here. As 
it has been shown in some SUSY realisations of flat or 
warped 5D models [l3|, at the massive level (i.e. not for 
the zero mode) the typical amount of SUSY which is left 
after the compactiflcation to four dimensions is higher 
then J\f = 1. There is the possibility that more than one 
tower of KK gravitinos is present, and the interaction 
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FIG. 1: Limits on the 5D Planck mass plotted against the 
coupling constant g^. The solid line is eq. II25I I. above it there 
is no -ff tx p regime after inflation, and no gravitational particle 
production. The lower dashed line is the constraint coming from 
the lowest KK modes: the region below the line is excluded (see 
eq. I|22| l and its relative counterpart accounting for a low reheat- 
ing temperature, eq. ((27J). The upper dashed line corresponds to 
eqs. ||23|I and l|26| l. obtained by demanding sufficiently low entropy 
release during BBN. The region below the dot-dashed line is ex- 
cluded instead thanks to the limits l|24l I25| l. 



and the axion decay constant /4, is 



1 



/4 



(29) 



Here = 92 MeV and to^ = 135 MeV are respectively 
the pion decay constant and mass, while z = mu/rrid « 
0.6 is the up to down quark mass ratio. 

The 4D theory is able to provide certain bounds on the 
axion decay constant using cosmological and astrophys- 
ical arguments, which depend on the parameters of the 
exact axion model, but the allowed range is essentially 
(see e.g. ref. [l9j ) 



10** GeV < /4 < 10'^ GeV 



(30) 



These limits imply a very light axion field, much lighter 
than 1 eV, as can be seen directly by use of which 
makes its production through gravity irrelevant. 

Once more, the scenario is radically different if the 
axion field lived in 5D. A specific 5D axion model was 
pushed forward in [20| . its warped counterpart is found 
in (2]| . In their model the anomalous interaction term, 
after dimensional reduction, reads 



strengths of gravitinos belonging to different towers can 
be different. That in turn means that, as shown in Q, 
(7* ~ g^,{MSSM) + (93/2 + 32 )r/ Am. which is even more 
favourable for gravitational particle production. Indeed 
higher energy density corresponding to the same tem- 
perature brings down the freezing out temperature (or 
can produce higher masses for the same temperature, as 
the Universe expands faster), strengthening the results 
outlined so far. 



KK Axions 

Axions are one more example of a weakly interact- 
ing particle that may endanger the successful contact 
between the standard cosmological model and observa- 
tional data. In the standard 4D theory, the axion pseu- 
doscalar field is introduced in order to solve the strong 
CP problem [l^, that is, to dynamically adjust the QCD 
O parameter to zero as a result of spontaneous symmetry 
breaking. 

The anomalous coupling of the axion to gluons, which 
breaks the Peccei-Quinn (PQ) symmetry is 



9 



32^2 



(28) 



where g is the strong coupling constant, the 4D de- 
cay constant, = G'^'^G'^^, with G^^ the gluonic field 

strength tensor and G^^ its dual, ^ a model dependent 
dimensionless quantity, and a is the actual axion field. 
The relation between the zero temperature axion mass. 
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•327r2 



/4 



(31) 



where ro = 1 and r„^o — V^, and a„ are the axion KK 
modes. These states are not the true mass eigenstates 
(see the discussion in [1^), which are needed for the scope 
pursued in this paper, but it can be easily shown that 
for high n these a„ come closer and closer to the actual 
mass KK eigenstates. Hence, each KK axion excitation 
couples to ordinary matter (gluons as well as photons) 
through the same 4D decay constant, whereas the 5D 
original axion field couples at a 5D PQ breaking scale /s 



r 
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32^2 



G2 



(32) 



where the two PQ scales are related by = 2nMr^f^ in 
the case of just one extra dimension. 

It is useful to divide the KK axion tower into two 
bands, the long-living modes and the unstable ones. The 
latter include axions which could decay during and after 
BBN. However, unlike the gravitino, in this case no con- 
straints can be derived for higher masses (lifetime shorter 
than 10~2 seconds) just by looking at BBN and DM den- 
sity, since there is no symmetry that protects the axion 
from decaying into photons or gluons, and there will con- 
sequently be no dangerous relics as long as the decay 
happens before the BBN. The decay rate will depend ex- 
plicitly on the mass of the axion and on /4. In order to 
make the analysis traceable, the 4D PQ scale will be fixed 
to /4 ~ 10^° GeV. In this work we are interested in grav- 
itational creation of particles, while other mechanisms as 
well as limits obtainable from astrophysical arguments 
are discussed in 1231. 
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The first band is that for which the axions, ignoring 
possible transitions between different KK modes, are sta- 
ble on the Universe's age timescale. The (n-th) axion 
decay rate a —s- 27 is 

^ 256^^ if' ^^^^ 

where a is the fine structure constant and C is a model 
dependent numerical coefhcient close to 1. Since the 
channel has a branching ratio of order 1, for C = 1 le 
lifetime is 

r = 2.10^°fg^Vf / Vyr (34) 
\ ma ) \ IQi" GeV ) ^ ^ ' 

and KK states with masses lower than about 25 keV have 
a lifetime exceeding the age of the Universe. The total 
number of these states is TV ~ 2 • lO^^/Mg. In com- 
plete analogy with the gravitinos, their contribution to 
the density fraction has to be set below the DM one, 
which means that 

^Ji^ - 4 • lO^'A/g"^^/'* < 0.12, 
A/5 > 11 TcV. (35) 

This result is not very stringent because, despite their 
weak interactions, only the very light axions are able to 
survive for so long, which translates in a poorly relevant 
gravitationally produced relic density. This is especially 
true when compared to the bounds extracted from grav- 
itinos, for they (likely) possess 7?.-symmetry, whose ana- 
logue is absent in the 5D axion scenario. 

A more interesting bound is that derived from the sec- 
ond band, which includes axions whose decay would al- 
ter the primordial abundances of light elements. Having 
fixed the 4D PQ scale, the masses associated to a given 
lifetime are automatically fixed as well ([33]) . 

Now, recalling that /4 ~ 10^ TcV, the states corre- 
sponding to a lifetime of about a second have masses 
around 20 GeV. These axions are strongly constraining 
if the decay branching ratio into hadrons is close to 1 [2^ , 
in which case they would efficiently alter the neutron to 
proton ratio during the early stages of BBN, affecting the 
final *He abundance. These qualitative considerations, 
once put into numbers, lead to M5 > 1.8 • 10^ TcV. 

If the hadronic channel is instead suppressed, then the 
main effect is the photodissociation of light elements at 
times t ~ 10^ s by decay of axions with mass about 
100 MeV. The aforementioned potential problem with 
such axions can be avoided if the 5D Planck mass is 
pushed beyond 2 • 10^ TeV. 

Notice that in both cases the corresponding mass gaps 
are tiny on the scales of interest (0.1 to 20 GeV), thereby 
justifying the approximations adopted in the preceding 
discussion. 

This is the most interesting limit coming from gravi- 
tational particle production of a 5D axion alone. Lastly, 
since the 4D PQ scale has been fixed, this lower limit on 



M5 leads to a lower limit on the 5D axion decay constant, 
namely /g > 10* GeV. 

One more comment is in order here. Since KK axions 
have weak but not gravitationally weak interactions, it is 
natural to ask whether they are in thermal equilibrium 
and when, for if that were the case then the gravitational 
abundances will be triggered to the usual thermal ones, 
erasing any information about their original production 
epoch. Consider for instance the 20 GeV KK excitation 
and the dominant a ^ 27 process. Roughly speaking, by 
comparing the rates of decay and of expansion H it 
is possible to obtain an approximate freezing out temper- 
ature at which this particular interaction becomes highly 
improbable. In the case under investigation one sees that 
when M5 > 1.3 • 10^^ TeV, then thermal equilibrium is 
established for Tt > Tg, > 1.4 • lO-^^Ml, i.e. during the 
H ^ p regime. Were the 5D Planck mass smaller, the 
inequality T^'y ^ H would return (this time in standard 
expansion regime) T^q > 1.7- 10^ TeV, but this tempera- 
ture is higher than Tt calculated for A/g < 1.3 • 10" TcV: 
there is never thermal equilibrium below the transition 
temperature. The constraints derived above are there- 
fore valid, the freezing out temperature for gravitational 
interaction being smaller than that for which thermal 
equilibrium would be realised. 



IV. CONCLUSION 

In this paper we have considered the phenomenology 
of gravitational particle production in the context of 
braneworld cosmology, where an epoch of non standard 
expansion [H ~ p) renders the mechanism appetible. 
Based on the results of 0, we have extended the anal- 
ysis to bulk fields that appear to a 4D observer as KK 
towers. As expected, the potential dangerousness of the 
produced abundances is able to strongly constrain the 
parameters of the extra dimensional model, confirming 
the well known fact that cosmology tends to dislike very 
low (e.g. TeV) gravity scales. 

Since the crucial ingredient which makes gravitational 
particle production a cosmologically important mecha- 
nism is a period of braneworld expansion, this mecha- 
nism is intrinsically entwined with the early Universe's 
inflationary era, for the temperature (or the Hubble pa- 
rameter) at the end of the reheating process depends, 
although the details are not fully understood, on the dy- 
namics of the infiaton field. 

Once a specific model of inflation is chosen, and a re- 
heating scale is computed, the abundances of gravitation- 
ally produced particles are automatically obtained, there- 
fore providing some information on the parameters of in- 
flation. This is so because there is a whole KK tower of 
particles accessible by gravitational production, and the 
reheating temperature resulting from inflation dictates 
the upper cutoff for the actual available modes. Hence, 
the results hitherto described have interesting implica- 
tions in braneworld inflation model building. 
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Two particular cases have been analysed in details, 
gravitinos and axions. This specialisation is required 
since, even though the production mechanism is basi- 
cally insensitive on the properties of the generated field 
(the Equivalence Principle), it is their interactions with 
SM fields that is essential in determining their late time 
history, and consequently the limits that can be inferred 
from it. 

If the 5D model is supersymmetrised then the presence 
of a KK tower (or more than one) of gravitinos strongly 
constrains the 5D Planck mass, at the level of 10^° TeV. 
This is especially true when the reheating process is effi- 
cient, allowing for extremely fastly expanding radiation 
dominated FRW Universe, which is the perfect environ- 
ment for gravitational particle production. 

A 5D axion is also able to put interesting bounds on 
M5 because, even though its mass (the gravity charge) 
is typically very small, and consequently not efficiently 
produced by gravitational interactions, the existence of 
KK axions of any mass makes it potentially relevant. The 
bounds obtainable by demanding that the gravitationally 
produced axions do not release too much entropy during 
the process of BBN, are around M5 > 10^ TeV (slightly 
model dependent) which, although not as strong as the 
KK gravitino bounds, is a fine example of how TeV scale 
gravity in the minimal extension of the SM employing a 
5D axion field can be ruled out by simple cosmological 
arguments. 

In general, although these explicit limits have been 
worked out only for two specific models, a similar anal- 
ysis can be performed for any beyond the SM bulk field 
with little extra effort, as long as these fields do not sub- 
stantially modify the brane Friedman equation [1] 

Concluding, gravitational particle production has a 
high phenomenological impact on braneworld cosmol- 
ogy, whose study provides yet another independent way 
of constraining the parameters of the extra dimensional 
spacetime. 
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APPENDIX A: MODELS OF BRANEWORLD 
INFLATION 

In this brief appendix the reasons why the focus has 
been chosen to be on chaotic inflation only are outlined. 
The main motivations are simply summarisable as, first, 
it is the one that (in the context of braneworld inflation) 



fits the WMAP constraints the best, and secondly, has 
the widest parameters space for which after reheating a 
period of brane expansion (H cx p) is likely to have taken 
place. 



1. Quartic Inflation 

The potential is V = Xcfi'^/A. The self coupling is con- 
strained by inflation to be A ~ 8 • 10~^^ 0, US- A possi- 
ble mass term for the inflaton field must be subdominant 
with respect to the quartic term for the duration of infla- 
tion, in order for this to be pure quartic inflation (and for 
the limit just reported to apply), and then this implies 
that the decay rate (using the same parametrisation as in 
the body of the paper) is going to be much smaller than 
the standard chaotic one, thereby implying further less 
restrictive limits (on the mass of the field in this case) . It 
is also likely that reheating completes during a standard 
p dominated regime. 



2. Exponential Inflation 

The potential is V ^ Vq cxp''°">'/"^ , where the initial 
vacuum energy is encoded in Vq. This model has several 
problems, first of all in braneworld scenarios it needs to 
be supplemented with a curvaton field in order to obtain 
the correct spectrum of the perturbations (in particular 
the scalar-to-tensor ratio) Secondly it is not inter- 

esting as far as gravitational particle production is con- 
cerned because its reheating (that takes place through 
gravitational decay of the inflaton) is very inefficient and 
almost always occurs in standard cosmological expansion 
era [H]. 



3. Natural Inflation 

The potential is ^ = [1 -I- cos{(f>/ f)], where / is the 
scale associated with the pseudo Nambu Goldstone bo- 
son and VA is the scale of inflation (vacuum energy) . 
One can obtain successful inflation on a brane with this 
potential provided that the following requirements are 
fulfilled A ~ 0.1M| and / ~ 3 • lO^A/g. The mass 
of the inflaton field is given by ~ A// ~ 3 • 10~^M5 
and is slightly heavier that the chaotic scenario. How- 
ever in this case reheating should proceed through loop 
interactions whose rate is F ~ m^/f^ ~ 3 • 10~^^M5. 
This means that everything goes as in the chaotic case, 
only with fixed coupling constant ~ 10""'^^, but this 
in turn then means that both flxh^ and Yx are both at 
least reasonably small (see for instance eq. [25|l . hence, 
harmless. 
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4. F-term and D-term Hybrid Inflation 

These models have been shown to be not rcahsable 
in the context of brane inflation, due to the (in)famous 
77 problem [23|i that is, it is not possible to have a po- 
tential that is flat enough to let the Universe inflate and 
expand as long as needed: inflation is very short lived (in 
fact, this is a general problem in the context of supersym- 
metric realisations of inflation, in particular arising from 
string theory [1^). 



APPENDIX B: ANOMALOUS PARTICLE 
CREATION 

Conformally invariant fields which live in a conformally 
flat spacetime (i.e. a spacetime conformal to Minkowski, 
like all the 2D spacctimes or the FRW Universes) are not 
gravitationally produced: indeed, one can always per- 
form a conformal transformation and obtain field equa- 
tions which look like those in Minkowski spacetime, al- 
lowing to define unambiguously a vacuum state. Exam- 
ples of conformally invariant fields are massless fermions, 
massless and conformally coupled scalars and, but only 
in 4D, the electromagnetic field. The presence of a 
mass term breaks conformal invariance and is respon- 
sible for particle production. Nevertheless, this is rig- 
orously true only at the classical level, since quantum 
effects may change this conclusion, as consequence of the 
trace anomaly. 

Indeed, after renormalization, the trace of the expecta- 
tion value of the energy-momentum tensor still vanishes 
for free massless fields in flat spacetime, but in more gen- 
eral cases it is not so. The trace anomaly has been calcu- 
lated for self-interacting fields in Minkowski spacetime: 
a well known example is the trace of a gauge field [l3 



2g 



(Bl) 



where (3 = ^dg/dfi is the /3-function of the rcnormal- 
ized coupfing g{fi), while C^^ is the strength tensor of 
the gauge field. Such a result can be easily understood 
by noticing that loop diagrams of (T^) diverge and the 
gauge invariant regularization is not conformally invari- 
ant. The only other case in which the trace anomaly has 
been computed rigorously is that of a free field in curved 
spacctimes. In 4D, one finds [ll| 



p. / Ten 



I 



28807r2 



1 



-02 ( RC.PR''^ - gi?' 



a3aR + a^R-' 



where a,; are numerical coefficients which depend on the 
spin of the field. On the other hand, there is no trace 
anomaly for the case of an odd number of spacetime di- 
mensions, which means that conformally invariant free 
fields propagating in a conformally fiat 5D bulk would 
be still traceless at the quantum level and cannot be pro- 
duced by an expanding FRW Universe. 

As discussed in rcf. the amplitude of particle pro- 
duction in a conformally fiat gravitational field is propor- 
tional to the trace of the energy-momentum tensor and 
therefore the trace anomaly can make the production of 
massless particles possible even in conformally fiat space- 
times. If we consider the case of an SU (N) gauge field, 
we find that for any power law expansion of the scale fac- 
tor, that is a{t) oc f^, the rate of particle production per 
unit time and volume is roughly 0^{g)/t'^. We can there- 
fore expect that, in the very early Universe, gravitational 
particle production via trace anomaly could be more ef- 
ficient than the one due to the mass term. For example, 
if we indicate with tji the time of the Universe after in- 
flation, the energy density of massless bosons produced 
by the expanding Universe to the cosmological energy 
density ratio in braneworld cosmology is 



P/Pc 



(B3) 



which can be compared with the counterpart in standard 
cosmology [T2j 



pIpc 



(B4) 



For M4 ^ Afs the mechanism may be efficient. Neverthe- 
less, interesting phenomenological consequences are un- 
likely: the process produces massless gauge bosons and, 
including radiative corrections, charged particles. On the 
other hand, interesting signatures would require weakly 
interacting particles which are never in thermal equilib- 
rium after being produced gravitationally and whose de- 
cay during or after the BBN may leave observable effects. 
From this point of view, there are currently no theoreti- 
cally well motivated candidates which can do the job. 



(B2) 



The non-zero trace in eq. (|B2[) may instead be respon- 
sible for the gravitational production of sterile or very 
weakling interacting particles. However, it is typically 
a very small quantity and, even if at high temperature 
may be more important than the production via the mass 
term (e.g. for superlight particles), it is not likely that it 
can lead to interesting and observable effects. 
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